Zero Temperature Coefficient Quartz Crystals for Very 
High Temperatures 

By W. P. MASON 

(Manuscript Received Nov. 15, 1950) 

In order to determine the angles of cuts for low temperature coefficient crystals, 
the elastic constants of quartz have been evaluated in the temperature range from 
— 100°C to -f200°C. This has been done by measuring a series of rotated Y-cut 
crystals in the thickness shear mode and a series of rotated X-cuts in the longi- 
tudinal length mode. From the measurements, low temperature coefficients AT, 
BT, CT, and DT type crystals can be determined which have their temperature of 
zero temperature coefficient at any prescribed temperature. Calculations are given 
for the properties of crystals to operate at 200°C. The characteristics of an AT 
type crystal have been investigated experimentally, and the measured results are 
in reasonable agreement with the calculations. It is shown that there is a maxi- 
mum temperature of 190°C for which an AT type crystal can have a zero temper- 
ature coefficient. 

I. Introduction 

Most quartz crystals used to control the frequency of oscillators or time 
measuring devices are used in places where the ambient temperature does 
not exceed 60° to 70°C. The crystals are usually adjusted in angle so that 
they have a zero temperature coefficient at a temperature of about 80°C 
and they are temperature controlled at this temperature. However, a 
class of uses occurs for which the ambient temperature may be considerably 
higher and for these uses ordinary AT and BT crystals, for example, are 
not satisfactory. This is evident from Figs. 1 and 2 which show the fre- 
quency variations for these crystals over a temperature range from — 100°C 
to +200°C. For example, the flattest frequency temperature curve for the 
AT cut occurs at an angle of +35°18' rotation about the X axis from the 
Y cut. By going to +35°36' orientation about the X axis a minimum occurs 
at 100°C. For the BT cut shown by Fig. 2 the angle of -49°16' orienta- 
tion gives nearly a parabolic shape centered at 20°C. By changing the orien- 
tation to -47°22' the parabola centers at 75°C. 

Hence if one wishes to raise the temperature for which the zero tempera- 
ture coefficient occurs he has to increase the rotation about X for the AT 
cut and decrease it for the BT cut. The amount needed for either orienta- 
tion can best be determined by evaluating the elastic constants as a func- 
tion of orientation and temperature, and that is the main purpose of this 
paper. The results are applied to determining the best angles of orientation 
for the AT, BT, CT, and DT type crystals to obtain zero temperature coeffi- 
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Fig. 1 — Frequency temperature characteristics of AT type crystals. 
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Fig. 2— Frequency temperature characteristics of BT type crystals. 



368 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1951 

cients for any arbitrary temperature. These calculated values have been 
checked experimentally for the AT type crystal and the angles and proper- 
ties are approximated by the calculations. It is shown that there is a critical 
angle of +36°26' which results in the highest temperature of 190°C for 
which it is possible to obtain a zero temperature coefficient AT type crystal. 

II. Evaluation of the Elastic Constants as a Function 
of Temperature 

A simple method for taking account of the temperature terms is to ex- 
pand the frequencies for the known cuts in powers of the temperature 
around some reference temperature. Since the data of Figs. 1 and 2 run 
from — 100°C to +200°C, a convenient temperature is 50°C. Then 

f = f m [l + ai (T - To) + a 2 (T - To) 2 + a 3 (T - T ) 3 + • • •] (1) 

Over this temperature range the frequencies measured can be accurately 
represented by terms including the cubic as the highest. If T is taken as 
50°C, equation (1) can be solved for the constants a\, o 2 , and a s and we 
find 

1 8 

— ■} 1/200° — /-100°J + ~ l/l25° — /-25°J 
fll = 



300 /bo- 
/-ioo° ~r /200 r 

(h = 



(2) 



a 3 = 



22,500 /«,. 

(/2OO /— 100°) 2[/i25° /— 25°J 

5,062,500/ 5 o° 



where the subscripts refer to the temperatures for which the frequencies 
are measured. If we apply these equations to the AT crystal cut at 35°18' 
and the BT at — 49°16', we find, for the frequencies, the equations 

/ AT = 1.661 X 10 5 [1 + .22 X 10- 6 (r - 50°) 

+ 8.9 X 10~ 9 (r - 50) 2 -f- 82 X 10~ 12 (r - 50) 3 + • • •] 

(3) 
/ BT = 2.547 X 10 5 [1 - 2.2 X 10- 6 (T - 50) 

- 55.5 X 10- 9 (r - 50) 2 - 73 X 10" 12 (r - 50) 3 + • • •] 

In order to obtain the frequency and the variation of frequency with 
angle, use is made of the equation for a thickness shear vibration 



f _ 1 . /cu _ _L A /c*6 cos 6 + eft sin 2 8 — 2cu sin 8 cos 8 

} 2ty -j uf p 



(4) 
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where / is the thickness of the crystal, p the density, the angle of the 
normal of the plate measured from the Y axis and cu , cu and c 6 s three of 
the seven elastic constants of quartz measured at constant electric field. 
This equation is valid for an infinite plate but is also a good approximation 
for a crystal whose cross-sectional dimensions are 30 to 40 times the thick- 
ness dimensions. Since there are three constants, the two measurements for 
the AT and BT cuts will give only two relations and we need a measurement 
for another angle. As discussed in Chapter X, Section 10.2 of "Piezoelectric 
Crystals and Their Application to Ultrasonics," ' the remaining cut can 
be obtained by measuring the thickness shear mode of a Y cut plate or the 
face shear mode of a Y cut crystal. The latter mode is considerably easier 
to dimension in order to obtain a frequency corresponding to the shear 
mode. Table I shows measurements for the frequency constant of a Y 
face shear mode for a crystal having the following dimensions: Length 
along the X axis is 36.86 mm, width along the Z axis = 7.625 mm; thickness 
along the Y axis is 0.990 mm. High harmonics were used and the frequency 
constant was obtained by dividing the frequency by the overtone order. 
This frequency is controlled by the cf\ elastic constant according to the equa- 
tion 

/ = J; J— or cl = 4/ m 2h 2 P (5) 

In calculating the cf 4 constant from the resonant frequencies measured, 
a correction is introduced by the temperature expansion constants of the 
crystal. This follows from equation (5) since h the frequency determining 
axis and p the density both change with temperature. From measurements 
quoted by Sosman 2 for the expansion along the Z axis and perpendicular 
to the Z axis, we find 

h = Zo[l + 7.8 X 10- 6 (r - 50) + 2.8 X 10" 9 (r - 50) 2 

- 1.5 x io- l2 (r - 50) 3 - •■•] 

(6) 

U, v = foil + 14.6 x io- 6 (r - so) + 6.3 x io- 9 (r - SO) 8 

- 1.9 X 10- 12 (r- 50) 3 + •••] 
Multiplying these together the volume expansion is 

v = Uyh = /o[i + 37 x io- 6 (r - 50) + 15.8 x io- 9 (r - 50) 2 

- 5 x io- l2 (r - 50) 3 + •••] 

1 Pizoelecteric Crystals and Their Application to Ultrasonics, W. P. Mason, D. van 
Nostrand Co. 1950, Section 10.2, page 204. 

2 Sosman, The Properties of Silica, Chemical Catalogue Co., 1927, pp. 386, 387. 
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Since the density is the inverse of the volume, the square of the frequency 
constant for a crystal whose dimension is measured at 50°C must be multi- 
plied by the factor 



Uyh 






in order to correct for the effect of temperature expansion on the elastic 
constant. This correction is shown by the third column of Table I. The 
fourth column is then the value of cf 4 for the various temperatures. The 
fifth column shows the values of the a x , a 2 and a 3 constants for the tem- 
perature variation of cf 4 ■ 

Table I evaluates one of the elastic constants of the frequency equation 
(4). To evaluate the other two constants, use is made of the frequency 

Table I 



Temperature °C 


Frequency Con- 
stant Kilocycle 
Centimeteis 


Correction for 
temperature ex- 
pansion 


c u dynes 

cm 2 
xl0->° 


E 

Constants for c 4 < equation 


-100 
-25 

50 
+125 

200 


237.22 
236.26 
235.07 
233.56 
231.78 


1.0029 

1.0016 

1.000 

0.9984 

0.9965 


59.82 
59.26 

58.58 
57.75 
56.78 


a, = - 171 X 10-« 
a 2 = - 212 X 10- 9 
a 3 65 X 10-" 



constants for the AT and BT cuts given by equation (3). Over a tempera- 
ture range the thickness / is given by the equation 

/ = /oS cos 2 4- Z° sin 2 0) (8) 

where /„ and l z are the values of unit lengths along the Y and Z axis ex- 
pressed as a function of temperature. Inserting the values of (6) and (7) 
in equation (3), the elastic shear constants for the AT and BT cuts become 

4* (AT) = 2.924 X 10"[1 - 12 X 10~ G {T - 50) 



4- 12.8 X 10- 9 (r - 50) 2 4- 172 X 10" 12 (r - 50) 3 + • • •] 
C66 K (BT) - 6.877 X 10 u [l - 20 X 10~ 6 (r - 50) 



(9) 



- 176 X 10- 9 (r - 50) 2 - 238 X 10" 12 (r - 50) 3 4- 

From equation (4), the frequency equation, we have 

Cb*(BT) = 0.4258c£ 6 + 0.5742 cf 4 4- 0.9890 cf 4 

cif (AT) = 0.6661 cf 6 + 0.3339 c? 4 - 0.934 cf 4 

Since cf 4 is already known, tke two equations can be solved for c't and 
Cu and we find 



(10) 
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& = 0.8892 c« (BT) + 0.9328 c£ (AT) - 0.8221 c^ 
cu = 0.6282 elf (BT) - 0.4016 c'£ (AT) - 0.2264 ef 4 

Inserting the values from Table I and equation (9) the three elastic con- 
stants become 

A = 58.58 X 10 ,0 [1 - 171 X 10~ 6 (r - 50) 

- 212 x io- 9 (r - 50) 2 - 65 x io- i2 (r - 50) 3 + • • •] 

A - 40.26 X 10 10 [1 + 168 X 10- 6 (r - 50) 

- 5 X 10- 9 (r - 50) 2 - 167 x io- i2 (r - 50) 3 + • • •] 

A = 18.20 x io'°[i + 90 x io- 6 (r - 50) 

- 270 X 10-°(T - 50) 2 - 630 X 10" 12 (r - 50) 3 -\ ] 

To determine the frequency and temperature coefficients for any angle, 
one substitutes the values of the elastic constants and the temperature ex- 
pansion coefficients in the frequency equation (4), which results in the 
expression 

P = lO lo [(3.8O2cos 2 4-5.526sin 2 - 3.426 sin cos 9) + 10~ 6 (r - 50) 
[668 cos 2 - 828 sin 2 - 336 sin cos - 13.5 sin 2 cos 2 0-46 

sin 3 cos 0] + 10- 9 (r - 50) 2 [395 cos 2 - 1160 sin 2 + 354 (13) 
sin 6 cos - 12 sin 2 cos 2 0-24 sin 3 cos 0] + 10~ I2 (r - SO) 8 
[310 cos 2 - 884 sin 2 0+1130 sin cos 0] + — ] 

III. Properties of AT and BT Cut Crystals Having Zero 
Temperature Coefficients at High Temperatures 
The process for obtaining high frequencies cuts of the AT and BT type 
that will have zero temperature coefficients at a high temperature — for 
example 200°C— is to substitute for T in equation (13) the value 

T = 200° + AT (14) 

Inserting this value in (13) and collecting the results in powers of AT, we 

find 

P X 10- 10 = [3.913 cos 2 + 5.373 sin 2 - 3.465 sin cos - 0.0023 
sin 2 cos 2 - 0.0075 sin 3 cos 9} + (AT) X lO" 6 
[807 cos 2 - 1236 sin 2 - 154 sin cos - 17 sin 2 9 
cos 2 - 53 sin 3 cos 0] + (AT) 2 X 10~ 9 [534 cos 2 9 (15) 
- 1558 sin 2 + 862 sin cos - 12 sin 2 cos 2 0-24 
sin 3 cos 9] + (AT)* X 10" l2 [310 cos 2 9 - 884 sin 2 
+ 1130 sin cos 0] 
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V. Evaluation for the Remainder of the Elastic 
Constants over a Wide Temperature Range 

In order to evaluate the remainder of the elastic constants measure- 
ments were made of the frequencies of a series of X-cut longitudinal crystals 
over the same temperature range from — 100°C to +200°C. 3 The longi- 
tudinal crystals measured had their lengths at —30°, 0°, +30° and +60° 
from the Y axis. For the four crystals measured the results are shown by 



10,005 




to,ooo_ 

70 80 90 100 110 120 130 140 150 160 170 180 190 200 210 220 230 
TEMPERATURE IN DEGREES CENTIGRADE 

Fig. 5— Frequency temperature curves for fundamental and third overtone for a crystal 
cut at 36°45' rotation. 

Table II. The analysis for f m = / 5 o° and the three constants d , o 2 and a 3 
are also shown. 

To correct for the temperature expansion coefficients, the increase along 
/* is given by the last equation of (16) while / ± 30° and / ± 60° are 

/ ± 30° = .25/, + .75/, = Z [l + 12.9 X 10~ 6 (r - 50) 

+ 5.42 X 10- 9 (r - 50) 2 - 1.8 X 10~ 12 (r - 50) 3 + 

/ ± 60° = .75/, + .25/* = / [1 + 9.5 X 10- 6 (r - 50) 

+ 3.68 X 10- 9 (T - 50) 2 - 1.6 x io- i2 (r - 50) 3 + 
Since the frequency of a long thin bar is given by the equation 

11 1 



(19) 



1 21 VpsX 



B> 

or 5o 2 = 



4/;/ 2 p 



(20) 



introducing the length correction from (19) and the density correction 
from (7) one can correct for the effect of temperature expansion. 

3 These measurements were made by T. G. Kinsley. 
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Applying these corrections to the frequency equations of Table II the 
resulting compliance constants become 

5 2 Vx) = A = 1-271 X 10- 12 [1 + 16.5 X 10- 6 (r - 50°) 

4- 58.5 X 10- 9 (r - 50°) 2 + 33 X W~ 12 (T - 50) 3 + ■ • •] 
5? 2 ( + 3o-x) - 0.8159 X 10-' 2 [1 + 96.4 X 10- 6 (r - 50) 

+ 276.5 X 10- 9 (r - 50) 2 + 219.4 X 10- ,2 (r - 50) 3 + • • • ] 
j&HMj = 1.402 X 10- 12 [1 + 114.4 X 10- 6 (r - 50) 

+ 178 X 10- 9 (r - 50) 2 - 91.6 X 10~ 12 (r - 50) 3 + • • •] 
f&ttra-z) - 0.8614 X 10- l2 [l + 186.4 X 10-«(r - 50) 

+ 302.2 X 10- 9 (r - 50) 2 + 385.3 X 10- 12 (r - 50) 3 + • • •] 
The equation for the compliance constant si' for an X-cut crystal at an 
angle 6 from the Y axis has been shown to be 

522 s = A cos 4 6 + s$ sin 4 d - 2su cos 3 d sin 

(22) 
+ (2sn + tu) sin 2 cos 2 

Solving for the constants in terms of the compliances for the four angles 
measured we find 

K K E S 2 2 (+30°) "" $22 (-30°) . E _ E | « J 

$11 = $22(0°X)5$14 = T~S t S 33 — $22(0°) T ^$22(+60°) 

4 E 2 B /„ E , E \ 10 B 2 E ( n\ 

— ^$22(30°) — T$22(-30°) J U$13 + $44,) = ~ ^"$22(0°) _ - $22(60°) KA^J 

, 28, . 26 E 

-V -q-$22(+30°) T "g $22(-30°) 

Hence adding the results we find 

sn = 1.271 X 10~ 12 [1 + 16.5 X 10" a (r - 50°) 

+ 58.5 X 10- 9 (r - 50) 2 +33X l0~ i2 (T - 50) 3 + • • •] 
si, = 0.971 X 10" I2 [1 4- 134.5 X 10- fl (r - 50) 

+ 144 x io- 9 (r - so) 2 + 570 x io~ l2 (r - 50) 3 + • 

s F u = -0.4506 X 10~ 12 [1 + 139.5 X 10~ 6 (r - 50) 
+ 40 x io- 9 (r - 50) 2 - 54 X 10~ 12 (r - 50) 3 4- • • •] 
(25 f, + sfO = 1-785 X 10~ 12 [1 + 300 X 10~ e (T - 50) 

+ 460 X 10- 9 (r - 50) 2 - 98 X io-> 2 (r - 50) 3 + • • •] 

4 See "Piezoelectric Crystals and Their Application to Ultrasonics," page 204, equation 
10.26. 
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All the compliance constants are now determined except sf 4 , sf 3 and 
5i 2 . From the relations for a crystal in the quartz class 5 

E E 

B _ ^66 B _ r,/ B B\ C44 /nC\ 

Sii — B~~^ a 5 *«• _ a*u — ^12; — -5—5 j- 2 UjJ 

C44C66 Cu C44 Cqq — C14 

the remaining constants can be obtained. Inserting the values of Cu , c^ 
and Cu from (12), we find 

*4*4 = 1.986 X 10- l2 [l + 201 X 10- 6 (r - 50°) 

+ 200 x io- 9 (r - so) 2 - 26 x io- l2 (r - 50) 3 + • • ■] 

5? 6 = 2.89 X 10" 12 [1 - 138 X 10- 6 (r - 50°) 

- is x io~ 9 (r - so) 2 + 3 x io- i2 (r - so) 3 h — 1 

( 26 ) 

sf 3 = -0.1005 X 10" 12 [1 - 678 X 10-«(r - 50°) 

- 2110 X 10~ 9 (r - 50) 2 + 610 x io- l2 (r - 50) 3 + • • •] 

su = -0.174 X 10- 12 [1 - 1270 X 10~ 6 (r - 50°) 

- 575 X 10- 9 (r - 50) 2 - 215 X 10~ 12 (r - 50) 3 + • • •] 

It is sometimes desirable to use the c values as a function of tempera- 
ture. The remaining values can be obtained from the relations valid for 
quartz 5 



R 

r, B S33 . S44 n B $88 S44 K Siz E SlX -\~ $12 /r t ' 7 \ 

zcu — — t — ; tew — — - ; C13 — — — ; c 33 = \ii) 

a p <x p a a 



a o a -a n d r. 

•S33 , S44 r>,E __ S33 544 _ _B SlS . E S\\ "f- 5l2 

where 



B f B . E \ r. E 2 a E 1 B B \ n B 2 

a = 533 (S11 + Su) — 2^i3 ; p = 544 (^ii — Su) — Isu 
c 3 B 3 = 104.8 X 10+ 10 [1 - 165 X 10~ 6 (r - 50) 

- 187 x io- 9 (r - 50) 2 - 410 X io- |2 (r - 50) 3 + ■ • • 

tf t = 9.6 X 10+ 10 [1 - 510 x io- 6 (r - 50) 

- 2000 X 10~ 9 (r - 50) 2 + 600 X 10- |2 (r - 50) 3 + • • • 
cfi = 86.75 X 10+ I0 [1 - 53.5 X 10~ 6 (r - 50) 

- 75 x io- 9 (r - 50) 2 - 15 x io- l2 (r - 50) 3 + • • • 

cy, = 6.15 X 10'°[1 - 3030 X 10- 6 (r - 50) 

- 1500 x io- 9 (r - 50) 2 + 1910 x io- i2 (r -50) 3 + • • 

B See Piezoelectric Crystals and Their Application to Ultrasonics, page 207. 



(28) 
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VI. Predicted Angles for CT and DT Face Shear 
Crystals 

The other two cuts of primary interest for frequency controlled oscillators 
are the CT and DT low frequency face shear modes. An exact solution 
for the frequency vibration of a face shear mode has not yet been obtained, 
but Hight and Willard 67 have pointed out an empirical relation that agrees 
with the measured frequencies over the entire range of angles of rotated Y 
cut crystals. This relation is for a square crystal 

/ = ip /JAp (29) 

where I is one edge dimension and sf B the shear elastic constant pertaining 
to the face shear mode. In terms of the orientation angle 6 

s^ = su cos 2 + su sin 2 + 4jf 4 sin cos (30) 

Introducing the values of sf 4 , ;?< and su from equations (24) and (26) 
the frequency becomes 

fi x i<r 10 = 

1A27 

[(1.986 cos 2 + 2.89 sin 2 - 1.802 sin cos 0) 

+ (399 cos 2 - 398 sin 2 - 251.5 sin cos 0) (31) 

X 10~ 6 (r - 50) + (397 cos 2 0-52 sin 2 
- 72 sin cos 0) X 10" 9 (r - 50) 2 4- (-52 cos 2 
4- 8.7 sin 2 + 98 sin cos 0) X 10" 12 (T - 50) 3 + • • • 

Since the formula is very approximate the small correction due to tem- 
perature expansions has been neglected. With this equation the indicated 
angles for zero temperature coefficient — which are obtained by setting the 

6 A Simplified Circuit for Frequency Substandards Employing a New Type of Low- 
Frequency Zero-Temperature-Coefficient Quartz Crystal, S. C. Hight and G. W. Willard, 
Proc. I.R.E., Vol. 25, No. 5, pp. 549-563, May 1937. The factor 1.23 agrees better with 
experiment than the value 1.25 given in the paper. 

7 Since this paper was written a much more nearly exact solution of a face shear mode 
vibration has been obtained by R. D. Mindlin and H. T. O'Neil. This solution is an ex- 
tension of the thickness shear vibration of a crystal published by Mindlin (Journal of 
Applied Physics, probably March issue 1951). For a square plate there are two solutions 
which are very close in frequency. For case A which corresponds to / of equation (29) ly- 
ing along the X direction the empirical factor F becomes 

F = 1.2718 - .03471g - .03727g* where g = « — , 

and sw' and s 55 ' are the elastic compliances corresponding to the rotated cuts. For the B 
case which corresponds to / lying along z' the same formula holds but g = 7^ —,• 
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multiplier of (T — 50) equal to zero and solving for the rotation angles 
0i and 02 — are 

0i = +36°20' and 2 = -53°50' (32) 

as compared to the experimental values of +38°20' and —52°, which repre- 
sents a shift of about +2° orientation for both angles. At these calculated 
angles the frequencies are within about 1.5 per cent of the experimental 
values and the curvature constants agree approximately with the measured 
values. 

To obtain the angles for any other temperatures, for example 200°C, 
we substitute 

T = 200 + AT (33) 

and obtain the expansion in powers of AT. For 200°C this results in 

fl x io- 10 = 

14^27 

[2.055 cos 2 + 2.829 sin 2 - 1.840 sin cos 0] 
+ [514 cos 2 - 413 sin 2 - 266 sin cos 0] X 10~ 6 Ar (34) 
+ [373 cos 2 0-48 sin 2 - 28 sin cos 0] X 10~ 9 (Ar) 2 
+ [-52 cos 2 + 8.7 sin 2 + 98 sin cos 0] 

X lO" 12 ^) 3 ) 

The zero temperature coefficient angles are obtained by setting the coeffi- 
cients of AT equal to zero giving 

514 cos 2 - 413 sin 2 - 266 sin 6 cos = (35) 

Solving for we find 

= +39°50 7 and -56° (36) 

If we add 2° to each of these in order to correct for the difference between 
the formula and the measured results at 50°C, the most probable angles for 
zero coefficients at 200°C are 

= +41°50' and -54° (37) 

At these angles the indicated frequencies and variations of frequencies 
with temperature should be 

= 41°51'; 

/ = — * 105 [1 - 63 X 10" 9 (AD 2 - 8 X 10- 12 (A7f + • ■ ■] 

= -54°; (38) 

/ = 2 ° 4 * 1Q6 [I - 14 X 10 ? (A7f + 8 X 10" 12 (Ar) 3 + • - •] 
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While these results are probably not very exact on account of the lack of 
an exact solution for the frequency of a face shear plate, they indicate the 
angles and approximate variations with temperature for high temperature 
plates. So far no experimental results have been obtained for crystals of 
this type. 
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